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Abstract 

As is well-known, the real quaternion division algebra EI is algebraically isomorphic to a 4-by-4 
real matrix algebra. But the real division octonion algebra O can not be algebraically isomorphic 
to any matrix algebras over the real number field R, because O is a non-associative algebra over R. 
However since O is an extension of H by the Cayley-Dickson process and is also finite-dimensional, 
some pseudo real matrix representations of octonions can still be introduced through real matrix 
representations of quaternions. In this paper we give a complete investigation to real matrix repre- 
sentations of octonions, and consider their various applications to octonions as well as matrices of 
octonions. 

AMS Mathematics Subject Classification: 15A33; 15A06; 15A24; 17A35 

Key Words: quaternions, octonions, matrix representations, linear equations, similarity, eigenvalues, 
Cayley-Hamilton theorem 

1. Introduction 

Let O be the octonion algebra over the real number field R. Then it is well known by the Cayley-Dickson 
process that any a S O can be written as 

a = a' + a"e, (1.1) 

where a', a" e H = { a = oq + aii + a2j + a^k | = fc^ = — 1, ijk = — 1, qq — 03 G M }, the real 
quaternion division algebra. The addition and multiplication for any a — a' + a" e, b = b' + b"e G O are 
defined by 

a + b=ia' + a"e ) + (&' + b"e) = {a' + b' ) + {a" + b" )e, (1.2) 

and 

ab^{a' + a"e){b' + b"e) = {a'b' - Wa") + {b" a' + a"W )e, (1.3) 

where b', b" denote the conjugates of the quaternions b' and b" . In that case, O is an eight-dimensional 
non-associative but alternative division algebra over its center field M, and the canonical basis of O is 

1, ei=i, 62 =j, 63 = fc, 64 = 6, es^ie, eg = je, 67 = fee. (1.4) 
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The multiplication rules for the basis of O are listed in the following matrix 
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(1.5) 



where Eg, = [1, 61, 



67]. Under Eq.(1.4) all elements of O take the form 
a = ao + ai6i H h 0767, 



(1.6) 



where oq — 07 € M, which can also simply be written as a = Re a + Im a, where Re a = oq. The conjugate 
of a is defined to be 

o = o' — a"6 = Reo — Ima. (1.7) 



This operation satisfies 



a + b = a + b, ah = ba 



(1.8) 



for all a, 6 G O. The norm of a is defined to be \a\ := Va^ = = y/aQ + al + h a^. Although 

is nonassociative, it is still an alternative, flexible, quadratic, composition and division algebra over 1 
that is, for all a, b G O, the following equalities hold: 



a{ab) 



-- a\ 
a 



{ba)a = ba? , {ab)a — a{ba) :— aba, 



a^ - 2(Re a)a+ lap = 0, (Ima)^ = -|Ima|' 
|a6| = |a||6|. 



(1.9) 
(1.10) 

(1.11) 
(1.12) 



As is well known, any finite-dimensional associative algebra over an arbitrary field F is algebraically 
isomorphic to a subalgebra of a total matrix algebra over the field. In other words, any clement in a 
finite-dimensional associative algebra over F has a faithful matrix representation over the field. For the 
real quaternion algebra H, it is well known that through the bijective map 



: a = ao + a\i + a^j + a-^k G 



is algebraically isomorphic to the matrix algebra 



M 



ao 
ai 

0,3 



ao 
as 
-a2 



-a2 
-as 

ao 
ai 



-a3 
a2 

-ai 
ao 



(1.13) 
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ao, ai, a2, as G 



(1.14) 



and (/)(a) is a faithful real matrix representation of a. Our consideration for matrix representations of 
octonions are based on Eqs.(l.l) — (1.3) and the result in Eq.(1.13). 
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We next present some basic results related to matrix representations of quaternions, which will be 
serve as a tool for our examination in the sequel. 

Lemma 1.1 JTsf. Let a — -\- aii + a2j + a^k G H 6e given, where qq — 03 G M. Then the diagonal 
matrix diag( a, a, a, a) satisfies the following unitary similarity factorization equality 



Q 



where the matrix Q has the following independent expression 



ao ~ai 


-02 


-as 


ai ao 


-as 


02 




ao 


-ai 


03 -02 


ai 


ao 



i j k 
1 k -j 

-k 1 i 



-k j 



1 



which is a unitary matrix over . 



Lemma 1.2||I|l. Let a, 6 £ H, and A G M. Then 

(a) a = b 0(a) = 

(b) <P{a + b) = (f>{a) + (f>{b), (t>{ab) = </)(a)(^(&), (t>{\a) = X(j){a), 0(1) 

(c) a — jE4(j){a)El, where E4 [1, i, j, k] and E\ [1, — i, — j, —kY". 

(d) 0(a) (a). 

(e) 0(a-i) =0-i(a), i/a^O. 

(f) det [0(a)] = |a|4. 

We can also introduce from Eq.(1.13) another real matrix representation of a as follows 



T(a) := K(f{a)K = 



where K — diag( 1, —1, —1, —1 ). Some basic operation properties on r(a) are 

r{a + b)^T{a)+T{b), T{ab) ^ T{b)T{a), T{a)^T'^{a), 
det [0(a)] = |a|^ <j){a-^) = <j>-^{a) if a 7^ 0. 
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(1.15) 



(1.16) 



(1.17) 



(1.18) 
(1.19) 



Combining the two real matrix representations of quaternions with their real vector representations, 
we have the following important result. 

Lemma 1.3. Let x — xo + xii + X2j + x^k G H. and denote lif = [xo, Xi, X2, x^ , called the vector 
representation of x. Then for all a, 6, a; G H, we have 



and the equality 



ai = (j){a)~x' , xi = T{b)lif , axh — (j){a)T{b)l^ = T{b)(j){a)l^ , 



0(a)r(6) = r(&)0(a) 



(1.20) 
(1.21) 
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always holds. 
Proof. Observe that 



it — 0(a;)a4 , it ~ T{x)a'^ , 
We find by Lemma 1.1 and Eq.(1.2) that 

and 



[1, 0, 0, 0]. 



ax)a^ ~ (l){a)<j){x)aj — (l){a)lt , xt = T{xb)aJ = T{b)T{x)a'^ = T{b)lt , 



axb = a{xb) = <j){a){xb) = (f){a)T{b)lt , axb = (ax)b = T{b){ax) = T(6)0(a)'af . 
These four equahties are exactly the results in Eqs.(1.20) and (1-21). □ 

Lemma 1.4[|l|] [0. Let a, b, x e O be given. Then 

(a) Re (ab) = Re (ba), Re {{ax)b) = Re {a{xb)). 

(b) {aba)x = a{b{ax)), x{aba) = {{xa)b)a. 

(c) {ab){xa) = a{bx)a, {bx){ab) — b{xa)b. 

(d) (a, b, x) — —(a, x, b) — (x, a, b), where (a, 6, a;) — {ab)x — a(bx). 

2. The real matrix representations of octonions 

Based on the results on the real matrix representation of quaternions, we now can introduce real matrix 
representation of octonions. 

Definition 2.1. Let a = a' + a"e E O, where a' = ao + aii + a2j + a-^k, a" = 04 + a^i + a^j + ajk E H. 
Then the 8x8 real matrix 

<j){a') -T{a")Ki 
(t){a")K^ T(a') 

is called the left matrix representation of a over M, where K/^ — diag(l, —1, —1, —1). Written in an 
explicit form, 

oo ^0-1 ^o-i —0,5 ^o,Q —aj 

ao -03 02 -05 



(2.1) 



uj{a) = 



ai 



02 03 ao -fli 



04 ay 
-ay a4 



03 -02 



ai 



ao —07 ag —05 



-ag 
a 5 
04 



a4 05 ag ay oq — ai 



-02 -a3 



as —04 



ay 



ai ao as -a2 



ae 



-ay — a4 a^ 02 



-a3 ao 



ai 



ay ag —as — a4 03 a2 — ai 



(2.2) 



Theorem 2.1. Lei x = xq + x\e\ 
vector representation of x. Then 



XyCy e 



id denote lit — [xq, xi, 



aS = Lu(a)l}f 



sy]"^, called the 



(2.3) 



/loWs /or a, a; e 
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Proof. Write a, x & O as a = a' + a"e, x = x' + x"e, where a', a", x' , x" € H. We know by Eq.(1.3) 
that ax = {a'x' - x^a") + {x"a' + a'V)e. Thus it follows by Eq.(1.20) that 

ai = 



x"o! + d'x' 



x"a! + o!'x' 

(j){a')^ - T{a")Kix" 
T{a')^' + (j){a")Ki~^ 

(j){a') -T{a")K4 
(t){a")Ki T(a') 



X 

7' 



as required for Eq.(2.3). 
Theorem 2.2. Let a e 



□ 



be given. Then 

aEs = Esu>ia), and E^a = u;{a)E^, 



(2.4) 



where '■— [I, ei, • • • , 67 ], and E^ :— [ 1, — ei, • • • , —67 ] . 

Proof. Follows from a direct verification. □ 

We can also introduce from Eq.(2.1) another matrix representation for an octonion as follows. 

Definition 2.2. Let a = a' + a"e = ao + aiei H + 0767 e O be given, where a', a" G H. Then we 

call the 8x8 real matrix 



u{a) := Ksu;^{a)K8 

the right matrix representation of a, where Kg = 
explicit form, 



(2.5) 



T(a') -(/-(a") 
0(a") T(^) 

diag(/r4, an orthogonal matrix. Written in an 
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(2.6) 



Theorem 2.3. Let a, x GO be given. Then 



= v{a)lt 



(2.7) 



holds. 



Proof. Write a, a; € O as a = a' + a"e, x = x' + x"e, where a', a" , x', x" G H. we know by (1.3) that 
xa = ( x'a' - aFx" ) + ( a"x' + x"ci/)e. Thus we find by Eq.(1.20) that 



xd = 



a"x' + x"a' 



a"x' + x"a' 



T(a')? - <p{a")x" 



T(a') -(j){a") 
<l>{a") T(^) 



as required for Eq.(2.7). 



□ 
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Theorem 2.4. Let a € O be given. Then 

aFs = Fsy^{a), and F^a = u'^{a)F^, (2.8) 

where Fs := [1, — ei, —67] and Fg := [1, ei, 67]-^. 
Proof. Follows from a direct verification. □ 

Observe from Eqs.(2.1) and (2.5) that the two real matrix representations of an octonion a = a' + a"e 

arc in fact constructed by the real matrix representations of two quaternions a' and a". Hence the 
operation properties for the two matrix representations of octonions can easily be established through 
the results in Lemmas 1.2 and 1.3. 

Theorem 2.5. Let a, b GO, A G K 6e given. Then 

(a) a = b oj{a) = ui{b). 

(b) ui{a + b) = u){a) + co{b), ui{Xa) = Xu){a), a;(l) = /g. 

(c) oj{a) = u)'^{a). 

Proof. Follows from a direct verification. □ 

Theorem 2.6. Let a, b € O, X €M. be given. Then 

(a) a = b <J=^ v{a) = iy{b). 

(b) i^{a + b) = i^{a)+i^{b), i^{Xa) = Xiy{a), i/(l) = Is. 

(c) i^{a) — v'^{a). 

Proof. Follows from a direct verification. □ 
Theorem 2.7. Let a G O be given. Then 

a=lE8u;{a)El and a = lFsiy'^{a)F;, (2.9) 
o o 

where E^, E^, F^ and F^ are as in Eqs.{2.4) and (2.8). 

Proof. Note that u}{a) and ^{a) are real matrices. Thus we get from Eqs.(2.4) and (2.8) that 

Es{E;a) = Es[uj{a)E;] = Esuj{a)E;, and F{F^a) = Fs[iy^{a)F^] = Fsiy^{a)F^. 
On the other hand, note that O is alternative. It follows that 

Es{E^a) = a — ei(eio) — • • • — 67(670) — a — e\a — ■ ■ ■ — e^a = 8a, 

and 

F^{F^a) = a — ei(eia) — • • • — 67(670) = a — e\a — • • • — 670 = 8a. 
Thus we have Eq.(2.9). □ 
Theorem 2.8. Let a GO be given. Then 

det [Lu{a)] = det [iy{a)] = \a\^. (2.10) 
Proof. Write a = a' + a"e. Then we easily find by Eqs.(1.21) and (2.5) that 



det [co{a)] = det [z^(a)] = 



T(a') -<?!)(a") 
Ha") t(^) 



= det [T{a')T{a') + 4>{a")4>{a") ] 



det[r(a'a') + (/<(a"a") 
det[|af /4 + |a'f A] 
{\a'\^ + \arf = \a\^. 
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as required for Eq.(2.10). □ 

Theorem 2.9. Let a & O be given. Then the two matrix representations of a satisfy the following three 
identities 

a;(a^) = a;^(a), i^(a^) = :^^(a), uj{a)i^{a) — v{a)uj{a). (2-11) 

Proof. Applying Eqs.(2.3) and (2.7) to the both sides of the three identities in Eq.(1.9) leads to 

Lo'^{a)'V = uj{a^)l)' , v'^{a)l) = v{a^)l) , i^{a)v{a)l) = v{a)uj{a)l} . 
Note that is an arbitrary 8x1 real vector when h runs over O. Thus Eq.(2.11) follows. □ 
Theorem 2.10. Let a G O be given with a ^ 0. Then 

uj{a~^) =ui~\a), and u{a~^) = u'^a). (2.12) 
Proof. Note from Eqs.(l.lO) and (1.11) that 

«"' = ^ = ^[2(Rea)-a] 

and 

- 2Rea+\af = 0. 

Applying Theorems 2.5 and 2.6, as well as the first two equalities in Eq.(2.11) to the both sides of the 
above two equalities, we obtain 

to{a-^) = -^[2{Rea)Is - io{a)], iy{a-^) = ^[2(Rea)/8 - i^{a)] 
\a\ \a\ 

and 

uj^{a) - 2(Rea)w(a) + \a\^Is = 0, i/^(a) - 2{Rea)i^{a) + [ap/g = 0. 

Contrasting them yields Eq.(2.12). □ 

Because O is non-associative, the operation properties uj{ab) — Lu{a)uj{b) and i>{ab) ~ v{b)i>{a) do 
not hold in general, otherwise O will be algebraically isomorphic to or algebraically anti-isomorphic to 
an associative matrix algebra over R, this is impossible. Nevertheless, some other kinds of identities on 
the two real matrix representations of octonions can still be established from the identities in Lemma 
1.4(a)-(d). 

Theorem 2.11. Let a, b G O be given. Then their matrix representations satisfy the following two 
identities 

Lu{aba) = Lu{a)uj{b)Lu{a), and v{aba) = v{a)i'{b)v{a). (2-13) 
Proof. Follows from applying Eqs.(2.3) and (2.7) to the Moufang identities in Lemma 1.4(b) . □ 
Theorem 2.12. Let a, b G O be given. Then their matrix representations satisfy the following identities 

oj{ab) + uj{ba) = uj{a)uj{b) + uj{b)uj{a), (2-14) 

v{ab) + v{ba) = v{a)v{b) + v{b)v{a), (2.15) 

uj{ab) + v{ab) = uj{a)uj{b) + v{b)v{a), (2.16) 

uj{a)v{b) + uj{b)v{a) = v{a)uj{b) + v{b)uj{a), (2.17) 

uj{ab) = uj{a)uj{b) + uj{a)iy(b) - iy(b)Lu{a), (2.18) 

i/(o6) = i/(&)i/(a) + oj{b)u{a) - u{a)u>{b). (2.19) 
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Proof. The identities in Lemma 1.4(d) can clearly be written as the following six identities 

{ab)x — a{bx) = —{ba)x + b{ax), {xa)b — x{ab) = —{xb)a + x{ba), 

{ab)x — a{bx) = —{bx)a + b{xa), {ab)x — a{bx) = —(xa)h + x{ab), 
{ab)x — a{bx) = —{ax)h + a{xb), {xa)b — x{ab) = —{ax)b + a{xb). 
Applying Eqs.(2.3) and (2.7) to the both sides of the above identities, we obtain 

[u{ab) — w{a)u){b) ]lt = [ —co{ba) + w{b)u{a) 
[ v{b)v{a) - v{ab) lot = [ -i/(a)z/(6) + u{ba) ]'^, 
[L'{b)u){a) — oj{a)v{b) ]Tt = [—i'{a)u){b) + u){b)i'{a) ]'^, 

[Lj{ab) — Lo{a)u{b) = [— j/(&)i/(a) + ^{ab) ] af , 
[uj{ab) — Lu{a)uj{b) ]lt = [—i>{b)Lu{a) + uj{a)v{b) 
[p{b)iy{a) - iy{ab)]lt = [-v{b)u{a) +w{a)p{b)]lt. 

Notice that is an arbitrary real 8x1 real matrix when x runs over O. Therefore Eqs.(2.14) — (2.19) 
follow. □ 

Theorem 2.13. Let a, b G O be given with a ^ 0, b ^ 0. Then their matrix representations satisfy the 
following two identities 

uj{ab) = i'{a)[uj{a)(j{b)]i'-'^{a), and u{ab) = oj{b)[v{b)v{a)](j~'^{b). (2.20) 

which imply that 

ui{ab) U!{a)uj{b), and v{ab) ^ u{b)v{a). (2.21) 

Proof. Applying Eqs.(2.3) and (2.7) to the both sides of the two identities in Lemma 1.4(c), we obtain 

uj{ab)y{a)lt = i'{a)uj{a)uj{b)lt , and v{ab)uj{b)l^ = uj{b)i'{b)i'{a)lt , 

which are obviously equvalent to Eq.(2.20). □ 

Note from Eqs.(2.3) and (2.7) that any linear equation of the form ax — xb = c over O can equivqlently 

be written as [w(a) — z^(6)]"5? = Ht, which is a linear equation over R. Thus it is necessary to consider 
the operation properties of the matrix Lu{a) — iy{b), especially the determinant of uj{a) — v{b) for any 
a, 6 G O. Here we only list the expression of the determinant of uj{a) — i/(&). Its proof is quite tedious 
and is, therefore, omitted here. 

Theorem 2.14. Let a, e O he given and define 5{a, b) := Lo{a) — v{b). Then 

det[5(a, b)] = |a-6|'^[s^ + (|Ima| - \lxnb\f][s^ + ( |Ima| + \l^b\f] (2.22) 
det [5{a, b)] = (s^ + |Ima + Imfep )2[s4 + 2s'^{\lma\^ + llmfeH + ( jlmap - |Im6|2 )\ (2.23) 

where s = Re a — Re 6. The characteristic polynomial of 5{a, b) is 

I A/8 - (5(a, 6) I 

= [{\-sf + \lma + Im6|2]2[ {X-sf + (|Ima| - \luvb\ f ][{\ - sf + { |Ima| + |Im6| f ]. (2.24) 
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In particular^ ifKea = Keb and |Inia| = |Im6|, but a ^ b, then 

rank(5(a, 6) = 6. (2.25) 

Theorem 2.15. Let a, 6 G O &e given with a ^ and b ^ 0. Then d{a, b) — uj{a) — v{b) is a real 
normal matrix over M., that is, d{a, b)d'^{a, b) ~ S'^{a, b)6{a, b). 

Proof. Follows from 

S{a, b) + S^{a, b) = uj{a) - ,y{b) + oj'^ {a) - ly'^ (b) 
— a;(a) — v{b) + uj{a) — v{b) 
= u}{a + a) - v{b + b) ^2{Rea-Reb)h. □ 

Theorem 2.16. Let a £ O be given with a ^ R. Then 

^^(a, a) = -4|Ima|^5(a, a), (2.26) 

and (5(a, a) has a generalized inverse as follows 

= «)• (2-27) 

Proof. Observe that S{a, a) — uj{a) — i'(a) = ijj(Inia) — i'(Inia) and (Ima)^ = — |Imap. Thus we 
find that 

S'^{a, a) = [u;(Ima) — :/(Ima) ]^ 

= [w^(Ima) — 2a;(Inia)i^(Ima) + i^^(Ima) ] 
= [a;((Ima)^) - 2uj{lma)iy{lma) + iy((Ima)^)] 
= -2[|Iniap/8 +cj(Ima)z/(Ima)], 

and 

S^{a, a) = — 2[ jlmap/g + a;(Ima)i^(Im a) ] [a;(Inia) — i^(Im a) ] 
= -4|Ima|^[w(Ima) - j/(Ima)] = -4|Im,a|^(5(a, a), 

as required for Eq.(2.26). □ 

3. Some linear equations over O 

The matrix expressions of octonions and their properties introduced in Section 2 enable us to easily 
deal with various problems related to octonions. One of the most fundamental topics on octonions is 
concerning solutions of various linear equations over O. In this section, we shall give a complete discussion 
for this problem. Our first result is concerning the linear equation ax — xb, which was examined by the 
author in 

Thoerem 3.1 iQ. Let a = ao + aici + • • • + 0767, b = bo + bici + • • • + 6767 £ O &e given. Then the 
linear eguation ax = xb has a nonzero solution if and only if 

Rea = Re6 and |Ima| = |Im6|. (3.1) 
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(a) In that case, ifbj^a,i. e., Ima + Im6 9^ 0, then the general solution of ax = xb can be expressed 

as 

X = (Im a)p + p{lm b), (3.1) 
where p e A{a, b), the subalgebra generated by a and b, is arbitrary, or equivalently 

X = Ai(Ima + Im6) + A2[|Ima| |Imfo| - (Ima)(Im6) ], (3.2) 

where Ai, A2 S K are arbitrary. 

(b) Ifb = a, then the general solution of ax = xb is 

X = xiei + X2e2 -\ h 0:767, (3.3) 

where xi xj satisfy aixi + 02X2 + • • • + 07X7 = 0. 

The correctness of this result can be directly verified by substitution. 

Based on the equation ax = xb, we can define the similarity of two octonions. Two octonions are said 
to be similar if there is a nonzero p G O such that a = pbp^^, which is written as a ^ 6. Theorem 3.1 
shows that two octonions are similar if and only if Re a = Re 6 and |Ima| = |Im6|. Thus the similarity 
defined here is also an equivalence relation on octonions. In addition, we have the following. 

Theorem 3.2. Let a, b G O be given with b^a. Then 

o ~ 6 -^=^ w(a) ~ u;(6). (3.4) 
Proof. Suppose first that a ^ b. Then it follows by Eq.(l.ll) that 

a^ - 2(Re a)a = -\af = -\bf = 6^ - 2{Reb)b. 
Applying Theorem 2.5(a) and Eq.(2.11) to the both sides of the above equality and we get 

Lo^{a) - 2{Rea)Lo{a) = Lo'^{b) - 2(Re6)w(6). 

Thus 

uj^{a) +U!{a)u;{b) - 2(Rea)a;(a) = u;'^{b) + u;{a)Lu{b) - 2(Re6)w(fo), 
which is equivalent to 

u){a)[u){a) + u){b) - 2(Rea)/8] = [u}{a)+u}{b) - 2(Re6)/8 ]a;(6), 

or simply 

u){a)u}{lma + Imb) = a;(Ima + lmb)co{b). 

Note that Ima + Im6 ^ 0. Thus a;(Ima + Im6) is invertible. The above equality shows that w{a) ui{b). 
Conversely, if ui{a) co{b), then tracea;(a) = tracea;(6) and \co{a)\ = |w(6)|, which are equivalent to 
Eq.(3.1). □ 

Next we consider some nonhomogeneous linear equations over O. 

Theorem 3.3. Let a, b G O be given with a ^ R. Then the linear equation ax — xa = b has a solution 
in O if and only if The equality ab = ba holds. In this case, the general solution of ax — xa = b is 

^ = .iT^ {ba-ab)+p- ^ (Ima)p(Inia), (3.5) 
4Imar Imor 
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where p E'Q is arbitrary. 

Proof. According to Eqs.(2.3) and (2.7), the equation ax — xa = b can equivalently be written as 

[Lj{a) - u{a) = 6{a, a)-^ = T. (3.6) 

This equation is solvable if and only if 5{a, a)6^{a, a) b = 5 . In that case, the general solution of 
Eq.(3.6) can be expressed as 

it = 5~ (a, a)l? + 2[ /§ — (5~ (a, a)6{a, a) J"^, 

where ^ is an arbitrary real vector. Substituting 

1 



5 (a, a) 



4|ImaP 



5{a, a), and 6 (a, a) = — 2[|Ima| + a;(Ima):^(Inia) ] 



in the; above two equalities and then returning them to octonion forms by Eqs.(2.3) and (2.7) produce 
the equality in Part (b) and Eq.(3.5). □ 

Theorem 3.4. Let a = ao + (XiCi H h aycy, b = bo + bie\ H h 6767 e O 6e given with a ^ R. Then 

the equation 

ax — xa = b 



has a solution if and only if there exist Aq, Ai G R such that 

b = Xo + Xia, 
in which case, the general solution of Eq.{3.7) is 



where xi — x-r satisfy 



Ai 

X = — + xici H h a;7e7. 



aixi H h arxr = — -Re 6. 



Proof. According to Eqs.(2.3) and (2.7), the equation (3.7) is equivalent to 



namely 



[u){a) — i'{a) = 6{ 


a, 


a)lt 


= b 




-2ai • • • -2a7 ' 




Xq 




" 60 ' 


2ai ••• 




X\ 




hi 


2a7 ••• 




. ^7 . 




. b7 _ 



(3.7) 
(3.8) 

(3.9) 
(3.10) 

(3.11) 



Obviously, this equation is solvable if and only if there is a Ai € K such that 

6i = Aiai, 62 = Aia2, 67 = Aia7, 

i. e., Imb = Ailma, which is equivalent to Eq.(3.8). In that case, the solution to xq is xq = and 
Xi — xy are determined by Eq.(3.9). □ 

Next we consider the linear equation 



ax — xo = c 



(3.12) 
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under the condition a ^ b. Clearly Eq.(3.12) is equivalent to 

[Lj{a) - u{b)]-^ = 6{a, h)lt = -$. (3.13) 

Under a ^ 6, we know by Theorem 3.3 that ax = xb has a nonzero solution. Hence 6{a, b) is singular 
under a b. In that case, Eq.(3.12) is solvable if and only if 

S{a, b)6^{a, b)"^ = -t, (3.14) 

and the general solution of Eq.(3.13) is 

-^ = 6-{a, b)-t + 2[Is-S-ia,b)S{a,b)]-^, (3.15) 

where ~jf is an arbitrary real vector. If a is not not similar to b. Clearly Eq.(3.13) has a unique solution 

lt = S-\a, b)-t (3.16) 

Eqs.(3.15) and (3.16) show that the solvability and solution of the octonion equation (3.12) can be 
completely determined by its real adjoint linear system of equations (3.13). Through the characteristic 
polynomial (2.24), one can also retern Eqs.(3.15) and (3.16) to octonion forms. But their expressions are 
quite tedious in form, and are omitted here. 

Another instinctive linear equation over O is 

a{xb) - {ax)b = c, (3.17) 

which is also equivalent to 

{ab)x — a{bx) = c, (3.18) 

as well as 

x{ab) - {xa)b = c, (3.19) 

because {ab)x — a{bx) = {ab)x — a{bx) = x{ab) — {xa)b hold for all a, b, x G O. Now applying Eqs.(2.3) 
and (2.7) to the both sides of Eq.(3.17), we obtain an equivalent equation 

[uj{a)i^{b) - v{b)Lo{a) \lt = (3.20) 

Here we set /i(a, b) = u){a)v{b) — i'{b)u}{a). Then it is easy to see that Eq.(3.20) is solvable if and only if 

/i(a, b)ii~{a, b)~$ = 

where /i~(a, b) is a generalized inverse of /x(a, b). In that case, the general solution of Eq.(3.20) is 

It = 11- {a, b)-t + [h- n-{a, b)iJ.{a, 6)]^, (3.21) 

where is an arbitrary real vector. Numerical computation for Eq.(3.21) can reveal some interesting 
facts on Eq.(3.17). The reader can try to find them. 

Theoretically speaking, any kind of two-sided linear equations or systems of linear equations over O 
can be equivalently transformed into systems of linear equations over R by the two equalities in Eqs.(2.3) 
and (2.7). Thus the problems related to linear equations over O now have a complete resolution. 
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4. Real adjoint matrices of octonion matrices 



In this section, we consider how to extend the work in Sections 2 and 3 to octonion matrices and use 
them to deal with various octonion matrix problems. Since octonion algebra is non-associative, the 

matrix operations in O is much different from what we are familiar with in an associative algebra. Even 
the simplest matrix miiltiplication rule A^A = AA^ docs not hold over O, that is to say. nniltiplication 
of matrices over O is completely not associative. Thus nearly all the known results and methods on 
matrices over associative algebras can hardly be extended to matrices over O. In that case, a unique 
method available to deal with matrices over O is to establish real matrix representations of octonion 
matrices, and then to transform matrix problems over O to various equivalent real matrix problems. 

Based on the two matrix representations of octonions shown in Eqs.(2.2) and (2.6), we now introduce 

two adjoints for a octonion matrix as follows. 

Definition 4.1. Let A = (ost) € O™^" be given . Then the left adjoint matrix of A is defined to be 



uj{A) = [oj{ast)] = 
the right adjoint matrix of A is defined to be 
u{A) = [u{ats)] = 



w(aii) 
zv(aii) 



Z^(aml) 



pSmxSn 



(4.1) 



pSnxSm 



and the adjoint vector of A is defined to be 
vec^ := [on'^, ■ 



>T — >T 

Ctml J ai2 , 



am2 



>T 



(4.2) 



(4.3) 



Definition 4.2. Let A = {Ast)mxn and B = {Bst)pxq are two block matrices over R, where Ast, Bst € 
]]^8x8 ^\^QTiy ^ijg igji g^jifj^ right block Kronecker products of A and B, denoted respectively by A^B and 
A^B, are defined to be 



A<SiB 



and 



AigiB = 



where 



AhQlB ■■ 

Ami QlB • • 
AQrBii •• 

A Qr Bpi ■ ■ 
AstBii 



Ast &lB = 



A Qr Bst = 



AstBpi 
AiiBst 

Ami Bst 



Aln ©L B 
Amn ©L B 

AQRBlg 

AQRBpg _ 

■ ■ AstBlq 

AstBpq 
AlnBst 

Bst 



x8nq 



(4.4) 



(4.5) 



X8q 



8mx8n 



(4.6) 



(4.7) 
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Noticing the equality (2.5), we see the two adjoint matrices co{A) and i^(^) of an octonion matrix A 
satisfy the foUowing equahty 

U{A) = Ksnto'^{A)K8m, (4.8) 



where 



Kst = diag(/r8, 



Ks), -fCg = diag(l, -1, -1), t = m,n. 



(4.9) 



It is easy to see from Eqs.(4.4) and (4.5) that the two kinds of block Kronecker products are actually 
constructed by replacing all elements in the standard Kronecker product of matrices with 8x8 matrices. 
Hence the operation properties on these two kinds of products arc much similar to those on the standard 
Kronecker product of matrices. We do not intend to list them here. 

We next present some operation properties on the two real matrix representations of octonion matrices. 

Theorem 4.1. Let A, B e O'"^", A G K 6e given. Then 

(a) A = B <J=^ uj{A) = uj{B) v{A) = v{B), i. e., uj and v are 1-1. 

(b) uj{A + B)=uj{A)+uj{B), and v{A + B) = v{A) + v{B). 

(c) u{\A) = \uj{A), and u{XA)=\v{A). 

(d) Lo{Im) = hm, and v{Im) = hm- 

(e) u){A*) = uj'^{A), and v{A*) = i''^{A), where A* = (o*^) is the conjugate transpose of A. 
Theorem 4.2. Let A e O™^" be given. Then 

1 



A = ^Esm^{A)Ej^, 



(4.10) 



where 



Est = diag( Es, ■■■ , Es), and Es = diag( 1, ei, 
Proof. Follows directly from Corollary 2.7. □ 



67 ), t = m, n. 



Since the multiplication of matrices over O is completely not associative, no identities on products of 
octonions matrices can be established over O in general. Consequently, no identities on products of the 
two kinds of real matrix representations of octonion matrices can be established. In spit of this, we can 
still apply Eqs.(4.1) and (4.2) to deal with various problems related to octonion matrices. Next are some 
results on the relationship of u){-), and vec(-) for matrices over O. 

Lemma 4.3. Let A e 0"^S B e O^^" and x G O be given. Then 

vec{Ax)=u}{A)lg' and vec{xB) = v{B'^)lt . (4.11) 
, a^Y and B = [bi, bnf. Then by Eqs.(2.3), (2.7) and Eqs.(4.1)— (4.3) 



Proof. Let A 

we find 



[ai, 



vec{Ax) 











uj(ai) 






uj{an)^ 




_ w(a„) 



X* = u){A)^, 



and 



vec{xB) = 



xbi 




' u{bi)lt ' 






xb^ 




_ I'ibn)'^ _ 




_ J^{bn) _ 



It = u{B^)lt. 



□ 
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Lemma 4.4. Let A e O™^", X e 0"^^ and a € O be given. Then 

yec{AX) = w{A)vecX and yec{Xa) = [i/(a)(g)/8n ]vecX = z/(a)0vecX. 

Proof. Let A = [Ai, ■ ■ A„] and X = [xi, ■ ■ Xn]'^- Then by Eq.(4.11) we find 

yec{AX) = vcc(^i.ti + |-A„x„) 

= vec(AiXi) + • • • + vec(A„x„) 
= uj{Ai)vecxi + ■ ■ ■ + Uj{An)vCCXn 

VCCXl 

= [oj{Ai), ■■■ , uj{An)] : =uj{A)vecX, 
as required for the first equality in (4.12). On the other hand, 



(4.12) 



vec(Xa) = 



xid 




i^(a)xi 


X^i^ 




via)x^ 



[i^(a)(8)/8n ]vecX = i'{a)(^vecX, 



as required for the second equality in (4.12). □ 
Lemma 4.5. Let B e O^^^ and X e 0"^p he given. Then 

vec(XS) = [!y(B'^)078„]vecX. 



(4.13) 



Proof. Let X = [Xi, ■ ■ ■ , Xp] and B = [bi, ■ ■ ■ , bp]^ . Then it follows from the second equality in 
(4.12) that 



Yec{XB) 



vcc(Xifei H h Xp&p) 

vec(Xifei) H h vcc(Xpfep) 

(i/(6i)(§/8„)vecXi H h (i^(6p)(§/8„)vccXj 

vecXi 

{[i^{bi), v{bp)]§>Isn) 



vecX„ 



iy{B^)§)Is„]vecX, 



as required for Eq.(4.13). □ 

Based on the above several lemmas, wc can find the following three general results. 
Theorem 4.6. Let A = (a^t) e O"^" and X e 0"^p be given. Then 

vec(^X) = [Isp^uj{A)]vecX. 

Proof. Let X = [Xi, • • • , Xp]. Then we find by Eq.(4.12) that 

vec{AX) = vcc[AXi, AXp] 

= [vec{AXi), • • • , vec{AXp) ] 

= [a'(j4)vecXi, uj{A)vecXp] 

= diag(w(A), uj{A))[\ecXx, vecXp] = [Isp(^co{A)]vecX, 



(4.14) 
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establishing Eq.(4.14). □ 

Theorem 4.7. Let B = (bst) e O*""' and X e ©"""^ be given. Then 

vec{XB) = [zy(B^)§/8„]vecX 

Proof. Let B = [Bi, • • • , Eg]. Then we find by Eq.(4.13) that 

vec(XB) = vec[XBi, XBg] 
vecXBi 



(4.15) 



vecXBa 



[iy{Bf)^Isn] yecX 



[.{Bf 



vecX = [i/(B^)g/8„]vecX, 



Proof. 



[z/(Bj)§/8„] vecX 
as rerquircd for Eq.(4.15). □ 

Theorem 4.8. Let A = (a^t) e O"^", B = (bst) G Op^«, and X e O"^^ 6e given. Then 
vec[{AX)B] = [p{B^)^u;{A)]vecX, and vec[A{XB)] = [uj{A)^p{B^)]vecX. 
According to Eqs.(4.14) and (4.15), we find that 
yec[{AX)B] = [u{B'^)§Ism]yec{AX) 

= [ v{B'^)%hm ] [ Isp%o^{A) ]vecX = [ v{B'^)%uj{A) ]vecX, 



i/(S^) ]vecX, 



(4.16) 



and 

Yec[A{XB)] = [l8p^w{A)]Yec{XB) 

= [l8p^u){A)][iy{B'^)^Isn]yecX = 

as required for Eq.(4.16). □ 

Theorem 4.9. Let A = (a^t) e O"^", X = {bst) e 0"^p, Y e 0«^" be given, and denote 
A^'^l *X = A{A ■ ■ ■ {AX) ■ ■ •)), and Y * A^''^ = ((••• {YA) ■ ■ ■)A)A. 

Then 

vec(A('=l *X) = [l8p^uj''{A) ]vecX, and vec{Y * A^''^) = [u''{A'^)^l8g ]vecY. (4.17) 

Just as the standard Kronecker products for matrices over any field, the three formulas in Eqs.(4.14) — 
(4.16) can directly be used for transforming any linear matrix equations over O into an ordinary linear 
system of equation over R. For example, 

[I§)uj{A) ]vecX = vecB, 
[i^(A^)(g)/]vecX = vecB, 
=^ [I^uj{A)lj{B) ]vccX = vecC, 
=> [iy{B^)iy{A^)®I]vccX = vecC, 
=^ [ iy{B'^)®Lj{A) ]vocX = vecC, 
=^ [uj{A)i^i^{B) ]vecX = vecC, 

» [I^uj{A) - iy{B)^I]vecX ^ vccC, 
B [iy{A^)^Lu{A) - co{A)(^iy{A^) ]vecX = vecB. 



AX = B ^ 
XA^B ^ 
A{BX) ^ 
{XA)B = C 
{AX)B = C ^ 
A{XB) = C ^ 
AX -XB=C ^ 
{AX) A - A{XA) 
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ThcorcticUy speaking, various problems related to linear matrix equations over the octonion algebra now 
have a complete resolution. 

Below are several simple results related to solutions of linear matrix equations over O. 

Definition 4.3. Let A e O"^"" be given. If its left adjoint matrix co{A) is invertible, then A is said to 

be completely invertible. 

Theorem 4.10. Let A = {agt) € O™^™ and B = (bst) G O™^" be given. If A is completely invertible, 
then the matrix equation 

AX = B, (4.18) 
has a unique solution over O. In that case, if the real characteristic polynomial of oj{A) is 

p{X) = A* + rt_iA*-i + • • • + riA + ro, (4.19) 
where ro is the determinant ofu}{A), then the unique solution o/£'g.(4.18) can be expressed as 

X = * B + rt-i(A(*-2| * 5) + • • • + r3A{AB) + r2AB + nB]. (4.20) 
Proof. According to Eq.(4.14), the matrix equation (4.18) is equivalent to 

[l8n^u}{A) ]vecX = vecB. (4.21) 

Because uj(A) is invertible, Igm^^i^) is also invertible. Hence the solution of Eq.(4.25) is unique and 
this solution is 

vecX = [l8n^oj{A)]~^vecB = [Ism^oj~'^{A)]vecB. 

Observe that 

w*(A) + rt-iw*-^(A) + • • • + riu){A) + rohm = 

holds. We then have 

1 



uj-\A)=-- [u;'-\A) + rt.iLo'-\A)+--- + r2io{A)+nIsm] ■ 



Thus 



l8n§)0J-^{A) = - — [/8„§W*-^(A) + rt_i(/8„Sa;*-2(A) ) + ••• + r2{l8n§)0j{A) ) + (ri/snS/sm ) ], 

and 

vecX= [Isn§'to(A)]-'^vecB = — (/8„0w*"^(A) )vecS + rt_i(/8„0a;*-2(A) )vecB 

_l 1_ r2{hn§>oj{A) )vecB + ri{hn®hm )veci?]. 

Retuning it to octonion matrix expression by Eq.(4.17), we obtain Eq.(4.24). □ 

Similarly we have the following. 

Theorem 4.11. Let A = {ast) & O™^™ and B = (bst) & O"^™ be given. If A is completely invertible, 
then the matrix equation XA = B has a unique solution over O. In that case, if the real characteristic 
polynomial of co{A) is 

p{X) = A' + rt-iA'-^ + • • • + riA + ro, (4.22) 
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then the unique solution of XB = A can be expressed as 

X = -^[B* + rt-i{B * + • • • + r3{BA)A + r2BA + nB]. (4.23) 

For simplicity, the two solutions in Eqs.(4.20) and (4.23) can also be written as 

X = L^^oB, X = BoR-^, (4.24) 

where L'^ and are, respectively, called the left and the right inverse operators of the completely 
invcrtiblc octonion matrix A. Some properties on these two inverse operators are listed below. 

Theorem 4.12. Let A e O™^™ be an completely invertible matrix, B e O'"^" and C G O"^™ be 
given. Then 

A{L-/oB)=B, A{L-/ o I^) = I^, (4.25) 

L-/ o {AB) =B, L-/oA = Im, (4.26) 

(Coi?-i)A = C, [Im o R^')A = Im, (4.27) 

(CA) oR^^ = C, Ao = (4.28) 

Proof. Follows from Theorems 4.10 and 4.11. □ 

We can also consider the inverses of octonion matrices in the usual sense. Let A e O™^™ be given. 
If there are X,Y£ Qmxm g^^]^ ^j^^^ xA = and AY = I„i, then X and Y are, respectively, called 
the left inverse and the right inverse of A, and denoted by Aj^^ := X and Aj^^ := Y. From Theorems 
4.10 and 4.11, we know that a square matrix of order m over O has a left inverse if and only if the 
equation [r^(A'^)(8)/8m ]vecX = veclm is solvable, and A has a right inverse if and only if the equation 
[l8m'i^^{A)]vecY = vecl„i is solvable. These two facts imply that the left and the right inverses of a 
square matrix may not be unique, even both of them exist. As two special cases, we have the following. 

Theorem 4.13. Let A G O™^™ be given. Then the left and the right inverses of A are unique if and 
only if A is completely invertible. In that case, if the real characteristic polynomial of uj{A) is 

p(A) = A* + rt_iA*-i + • • • + riA + ro, 

then the unique left and the unique right inverses A can be expressed as 

Al' = --[ + rt_a(*-2| + . . . + rsA{A') + r^A^ + n/^ ] , 
ro 

and 

A-^^ = --[ + ri_iAl*-2) + . . . + r^{A^)A + r-,A^ + ri/„ ], 

To 

where A^'\ := A{A{- ■ ■ {AA) ■ ■ •)) and A^"'^ := ((• • • {AA) ■ ■ ■)A)A. 

Proof. Follows directly from Theorems 4.10 and 4.11. □ 

Based on Theorems 4.10 and 4.12, as well as Eqs.(4.25) — (4.28), we can also derive the following two 
simple results. 

Corollary 4.14. If A € (□'"x™ is completely invertible, and ABi = ACi and B^A = C2A, then 
Bi — Ci and B2 = C2. In other words, the left and the right cancellation rules hold for completely 
invertible matrices. 
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Corollary 4.15. Suppose that A e O™^™, B e O"^" are completely invertible and C € O™^". Then 

(a) The matrix equation A{XB) = C has a unique solution X = {L^^ o C)R]^^. 

(b) The matrix equation {AX)B = C has a unique solution X = L~^(C o K^^, 
where and are the left and the right inverse operators of A and B respectively. 

Our next result is concerned with the extension of the Cayley-Hamilton theorem to octonion matrices, 
which could be regarded as one of most successful applications of matrix representations of octonions. 

Theorem 4.16. Let A e O"*^™ he given and suppose that the real characteristic polynomial of uj{A) is 

p{\) = A' + rt-iX'-^^ + • • • + nA + ro. 

Then A satisfies the following two identities 

+ ^ . . . ^ r:iA{AA) + r^A^ + nA + r^Im = 0, (4.29) 

AI*) + rt_iA'*-i) + • • • + r3{AA)A + ra^' + r^A + rolm = 0. (4.30) 
Proof. Observe that p[a;(A)] = 0. It follows that 

[ /8m0p[w(A)] ]vec/„ = 0. (4.31) 

On the other hand, it is east to see by Eq.(4.17) that 

vecA^^I = vec(A(''l * /„) = [Ism^io''{A) ]vecJ„, s = 1, 2, • • • . 

Thus we find that 

[l8m^p{uj{Aj)]vecIm 

= [Ism®oj\A) + r(_i(/8mSw*-i(>l) ) + ••• + ri{Ism^uj{A) ) + roihm^h-m ) ]vec J„ 

= (78m§w*(A) )vec/TO + rt_i(/8m§w*~^(^) )vec/m H h ri(/8mSw(A) )vec/m + roihm^hm )yeclm 

= vec A^*! + rt_ivec A^*"^! H h nvec A + rovec 

= vec[^(*l + rt_i^(*-i| + • • • + nA + rolm ]■ 

The combination of this equality with Eq.(4.31) results in Eq.(4.29). The identity in Eq.(3.30) can be 
established similarly. □ 

Finally wc present a result on real eigenvalues of Hermitian octonion matrices. 

Theorem 4.17. Suppose that A € O™^™ is Hermitian, that is, A* = A. Then A and its real adjoint 
Lu{A) have identical real eigenvalues. 

Proof. Since A = A*, we know by Theorem 4.1(e) that uj{A) = lu{A*) = a;^(A), that is, w{A) is a 
real symmetric matrix. In that case, all eigenvalues of uj{A) are real. Now suppose that 

Lo{A)X = X\, (4.32) 

where A e M and X e R^'^^^ Then there is unique Y e O™^^ such that vecF = X. In that case, it is 
easy to find by Theorem 4.1(a) and Eq.(4.12) that 

w{A)X = XX^ oj{A)vec Y = vec YX =^ vec (AY) = vec (FA) =^ AY = YX, (4.33) 
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which imphes that A is a real eigenvalue of A, and F is a eigenvector of A corresponding to this A. 
Conversely suppose that AY = YX, where A £ R, F G O™^^. Then taking vec operation on its both 
sides according to Eq.(4.12) yields 

u!{A)vecY — vec FA. 

This implies that A is also a real eigenvalue of uj{A) and vecY is a real eigenvector of U!{A) associcated 
with this A. □ 

The above result clearly shows that real eigenvalues and the corresponding eigenvectors of a Hermitian 
octonion matrix A can all be determined by its real adjoint Since uj(A) is a real symmetric 8m x 8m 

matrix, it has 8m eigenvalues and 8m corresponding orthogonal eigenvectors. 

Now a fundamental problem would naturally be asked: how many different real eigenvalues can a 

Hermitian octonion matrix A have at most? For a 2 x 2 Hermitian octonion matrix A ^ - , where 

a, c e M, its real adjoint is 

u:{A) = 

Clearly the characteristic polynomial of uj{A) is 

det( A/i6 - uj{A) ) = [ (A - a)(A - c) - \b\^ f. 

This shows that uj{A), and correspondingly A, has 2 eigenvalues, each of which has a multiplicity 8. 

The eigenvalue problem for 3x3 Hermitian octonion matrices was recently examined by Dray and 
Manogue |^ and Okubo They showed by algebraic methods that every 3x3 Hermitian octonion 
matrix has 24 real eigenvalues which are divided into 6 groups, each of them has multiplicity 4. Now 
according to Theorem 4.17, the real eigenvalues of any 3x3 Hermitian octonion matrix 



ah i^{b) 



A = 



ail ai2 ai3 

ai2 022 ^23 
ai3 023 033 



can be completely determined by its real adjoint 



;(A) = 



cj(oii) 
u;^(oi2) 

W^(oi3) 



Oil, 022, O33 G 



07(012) ^^(013) 
W(022) 0;(023) 
^^(023) ^(033) 



Obviously this matrix has 24 real eigenvalues and the 24 corresponding real orthogonal eigenvectors. 
Numerical computation shows that these 24 eigenvalues are divided into 6 groups, each of them has 
multiplicity 4, which is consistent with the fact revealed in 1^ and Moreover the 24 real orthogonal 
eigenvectors can also be converted to octonion expressions by (4.33). 

Furthermore, numerical computation reveals an interesting fact that the 32 real eigenvalues any 4x4 
Hermitian octonion matrix are divided into 16 groups, each of them has multiplicity 2; the 40 real 
eigenvalues of any 5x5 Hermitian octonion matrix are divided into 20 groups, each of them has multiplicity 
2. 

In general, we guess that for any m x m Hermitian octonion matrix with m > 3, its 8m real eigenvalues 
can be divided into 4m groups, each of them has multiplicity 2. 

As a subsequent work of Thereom 4.17, one might naturally ask hwo to establish a possible factor- 
ization for a Hermitian octonion matrix using its real eigenvalues and corresponding octonion orthogonal 
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eigenvectors, speak more precisely, for an m x m Hermitian octonion matrix A, how construct a complete 
invertible octonion matrix P (unitary?) and a real diagonal matrix D such that A = PDP~^ using 
its 8to real eigenvalues and 8m corresponding octonion orthogonal eigenvectors. However, this problem 
seems quite curious, because the number of different real eigenvalues of an Hermitian octonion matrix 
is more than its order. This problem is also quite challenging, because various traditional methods in 
associative matrix theory are not applicable to this non-associative case. 

As pointed out in Hermitian octonion matrices can also have non-real right eigenvalues. Theoret- 
ically speaking, the non-real eigenvalue problem of Hermitian octonion matrices may also be converted 
to a problem related to real representations of octonion matrices. In fact, suppose that AX — XX, where 
A e O and X e O™^^. Then according to Eq.(4.12), it is equivalent to 

uj{A)vecX = iy{X)§)vecX, 

or alternatively 

-diag(iy(A), iy{X)]vecX ^0. 
How to find iy{X) satisfying the equation remains to further study. 

Conclusions. In this paper, we have introduced two pseudo real matrix representations for octonions. 
Based on them we have made a complete investigation to their operation properties and have considered 
their various applications to octonions and matrices of octonions. However our work could only be 
regarded as a first step in the research of octonion matrix analysis and its applications. Numerous 
problems related to matrices of octonions remain to further examine, such as: 

(a) How to determine eigenvalues and eigenvectors of a square octonion matrix, not necessarily Hermi- 
tian, and what is the relationship of eigenvalues and eigenvectors of a octonion matrix and its real 
adjoint matrices? 

(b) Besides Eq.(4.29) and (4.30), how to establish some other identities for octonion matrices through 
their adjoint matrices? 

(c) How to establish similarity theory for octonion matrices, and how to determine the relationship 
between the similarity of octonions matrices and the similarity of their adjoint matrices? 

(d) How to consider various possible decompositions of octonion matrices, such as, LU decomposition, 
singular value decomposition and Schur decomposition? 

(e) How to characterize various particular octonion matrices, such as, idempotent matrices, nipoltent 
matrices, involutary matrices, unitary matrices, normal matrices, and so on? 

(f) How to define generalized inverses of octonion matrices when they are not completely invertible? 

and so on. As mentioned in the beginning of the section, matrix multiplication for octonion matrices is 
completely not associative. In that case, any further research to problems related matrices of octonions 
is extremely difficult, but is also quite challenging. Any advance in solving the problems mentioned 
above could lead to remarkable new development in the real octonion algebra and its applications in 
mathematical physics. 

Finally we should point out that the results obtained in the paper can use to establish pseudo matrix 
representations for real sedenions, as well as, in general, for elements in any 2"-dimensional real Cayley- 
Dickson algebras. 
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